The cat's cradle, stirring, and topological complexity 
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There are several physical situations in which the 'tangling' of a loop is relevant: 
the game of cat's cradle is a simple example, but a more important application 
involves the stirring of a fluid by rods. Here we discuss how elementary topology 
constrains the types of mappings that can occur on a surface, for example when the 
surface is the domain of a two-dimensional fluid. 
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In the children's game of cat's cradle, a loop of string is first anchored around both hands. 
The fingers of each hand are then inserted into the opposite hand's loop and pulled apart, 
resulting in the situation shown in Figure [T] The motion of the fingers allows us to predict 
in advance just how long a piece of string we will need: the more intricate the motion, the 
more string is required. In fact the length of string is not a bad first definition of 'topological 
complexity, ' and we will see how this simple concept can be applied to the study of dynamical 
systems, particularly in two dimensions. The game of cat's cradle also suggests a connection 
between topology and combinatorics: after all, the number of distinct cat's cradle shapes 
that can be made is essentially a difficult combinatorial problem. Throughout this short 
note, we will explore the interplay between topology and combinatorics, and see what it 
reveals about dynamical systems on surfaces. 

On the practical side of things, this 
sort of approach to topological complex- 
ity i s use ful when studying stirring de- 
v j ces l!]2i3] pjg Ure [2] shows a top-down view 

of such a device: the black disks are stir- 
ring rods immersed in a two-dimensional 
white fluid, assumed very viscous. The 
rods are moved following the prescription 
shown in the inset. The result is that an 
initial loop gets wrapped around the rods 
in a topologically complex manner, and is 

stretched exponentially (blue line). This is very good for the ultimate goal, which is to en- 
hance the chemical reaction rate or blending of solutes: the blue line represents the contact 
area between the two reagents. The reaction rate is limited by the growth of this contact 
area, so having it grow exponentially is desirable. 




FIG. 1: The cat's cradle. 



Two handles are better than one 

Figure [3] shows the double-torus S 2) beloved of low-dimensional topologists. The subscript 
2 on £2 stands for the number of handles, or genus. Though it isn't encountered very much 
in the real world, it is a surface that's complex enough to be interesting, unlike the humble 
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ordinary torus S\ which is a bit too trivial. We are interested in transformations of a surface 
S, i.e., continuous mappings of S to itself. Such a mapping, if also invertible, is called a 
homeomorphism of S. (We assume that the surface is orientable, and that the mapping 
preserves this orientation.) Now in topology one only cares about gross features, so we 
assume that two mappings are equivalent if they can be continuously deformed into each 
other: they are then called isotopic. An equivalence class of mappings under isotopy is called 
an isotopy class. The set of all such classes forms the mapping class group MCG(S') of S, 
with the group operation obtained from composition of mappings. 

The mapping class group of the torus Si is par- 
ticularly simple: after we allow for isotopy, all that 
matters is how many times a mapping <p : Si — ► Si 
has caused us to wind around each of the two pe- 
riodic directions of the torus. This winding is best 
illustrated by looking at the action of a mapping on 
two nontrivial disjoint closed curves drawn on the sur- 
face, as in Figure |4| The red loop around the thin 
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direction of the torus is labeled 



the blue loop 



FIG. 2: A stirring device with four rods. 
The inset shows the sequence of rod mo- 
tions, and the blue line is a blob of dye 
that was filamented by the rods. (From 
Ref. E) 



The mapping in Figure 



can be written as 



its linear action transforms the red loop to 



and the blue loop to 



which indicates the 



number of twists of each image loop around each pe- 
riodic direction of the torus. In fact, in general we 
have MCG(Si) ^ PSL 2 (Z) (the "projective" version 
of the group of invertible two-by-two matrices with 
determinant 1), since we can wind in positive as well as negative directions, and the action 
of 4> on these loops is commutative in the case of the torus. Note that since the loops initially 
only intersect once, then their images only intersects once as well. 

According to the celebrated Thurston-Nielsen 
(TN) classification theorem,^ there are only three 
types of isotopy classes: finite-order, reducible, and 
pseudo-Anosov. A mapping : S — > S is finite-order 
if <p k is isotopic to the identity map, for some inte- 
ger k > 0. The mapping is reducible if it leaves 
a set of disjoint closed curves invariant: the surface 
is carved up such that the pieces map to each other. 

Finite-order mappings are not interesting to us because they are too simple; reducible ones 
are not interesting because one can then reapply the TN theorem to the individual pieces 
(with boundaries) until the mapping is no longer reducible. 

It is the third type — pseudo-Anosov — that is most dynamically interesting. For such a 
mapping, in addition to leaving no nontrivial closed curve "untouched" (otherwise it would 
be reducible), the action on loops causes their length to grow exponentially.^ The asymptotic 
factor by which the winding numbers of loops is multiplied (under repeated iteration) is called 
the dilatation, A > 1. Its logarithm is the topological entropy. The topological entropy is 




FIG. 3: The double-torus, 52. 



FIG. 4: The action on loops of a mapping of the torus (their orientation is not shown). 



closely related to the growth rate of loops, as discussed for the cat's cradle and stirring device 
examples in the introduction. The positive entropy of a pseudo-Anosov mapping tells us 
that loops will grow exponentially under repeated iteration, with growth rate log A. 

For the torus Si, the TN classification is achieved merely by looking at the trace of 
an element of MCG(Si) ~ PSL 2 (Z). If the absolute value of the trace is greater than 2, 
then the mapping is pseudo-Anosov^ If less than 2, finite-order; If equal to 2, reducible.^ 
One can check this easily by using the Cayley-Hamilton theorem for two-by-two matrices 
of unit determinant. The mapping of the torus in Figure [4], for instance, is pseudo-Anosov 
with A ~ 2.618. 

Of polynomials and indices 

Now let us return to the double-torus 5*2. The mapping class group in this case is 
more complicated than for the torus, but we can get a handle on part of it by considering 
MCG + (S , 2) — the subset of the mapping class group that also leaves invariant a vector field 
(as opposed to a line field in the general case). In that case the mapping class group can be 
deduced from the action on the four loops drawn in Figure [3} so that MCG + (S , 2) C SL 4 (Z)P 
But this time we do not have equality: there are many matrices that do not correspond to 
an element of the mapping class group. But which ones? 

The matrices we want are in the symplectic group Sp 4 (Z)P One important feature of 
these matrices is that their characteristic polynomial 

P(x) = x 4 + C\ x 3 + C2 x 2 + C3 x + C4 

satisfies C3 = c±, c 4 = 1, so that the polynomial is reciprocal — the coefficients read the 
same from left to right as from right to left. The largest root (in magnitude) of P(x) is the 
dilatation A; for a pseudo-Anosov mapping it is always real. We appeal to the Lefschetz 
fixed point theorem, which in our case says 

L{4>) = 2 - Tr(0,) = Y, In d(<M- 

p6Fix(<£) 

Here, L(<f>) is an integer called the Lefschetz number of 0; 0* G Sp 4 (Z) is a matrix repre- 
senting the isotopy class of by looking at its action on the loops in Figure [3] Tx(0*) is the 
trace; Fix(0) is the set of isolated fixed points of 0; and lnd(0,p) is the topological index 
of at the isolated fixed point p. The topological index is given by the signed number of 
turns of a vector joining a point x and its image <fi(x) as x travels counterclockwise around a 
small circle enclosing the singularity (see Figure [5]). The index can only assume two simple 
values: it is 1 — N if fixes the 2N 'hyperbolic sectors' around a fixed point, and it is +1 
(independent of N) if 4> permutes the sectors.^ This is due to the severe restrictions imposed 
by continuity and orientation-preservation.^ 
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Now comes the first crucial ingredi- 
ent: the Cay ley-Hamilton theorem and 
Newton's formulas^ tell us how to find 
Tr(0*s) recurs i ve iy by looking at the poly- 
nomial P(x): 

fc-i 

Tr(^) = -^c m Tr(^ m )-A; Cfc , 

771=1 

where Ck = for k > n. Hence, we can cal- 
culate the Lefschetz numbers of all powers 
of directly from the polynomial. For ex- 
ample, for the polynomial 

P(x) = x 4 — x 3 — x 2 — x + 1 

with dilatation A ~ 1.72208, we find 



FIG. 5: A singularity with 12 hyperbolic sectors 
(N = 6). Each sector can map to itself (left, in- 
dex 1 — 6 = —5) or to one of two other sectors 
(right, index +1). The index is defined as the num- 
ber of turns of a vector joining x to <f>(x) as x travels 
counterclockwise around a small circle. 



L(0) = 1, L(0 2 ) = -1, L(0 3 ) = -5, L(0 4 ) = -5, L(0 5 ) = -14, L(0 6 ) = -25 . . . . 

The second crucial ingredient is that the index lnd(0,p) at a fixed point p is heavily 
constrained. Around a fixed point, the vector field that remains invariant under can have 
a singularity, near which the mapping's action can be divided into 2N hyperbolic sectors as 
in Figure [5j with A" even. Assuming the polynomial P(x) has a positive dominant root, then 
as mentioned above there are two cases: lnd(0,p) = 1 — A^<Oif0 fixes the sectors at p 
(left in Figure [5]); lnd(0,p) = 1 if it permutes the sectors (right in Figure [5]). For technical 
reasons,^ a permutation involves only N/2 of the sectors. 

Let's see what we can glean from this for the case of the polynomial above: on a surface 
of genus two, the Euler-Poincare formula dictates 

(2 - AW = 2X2 = -4 

p6Fix(</>) 

where \2 — ~ 2 is the Euler characteristic of the surface 5*2, and 2N p is the number of 
hyperbolic sectors at a fixed point p. 

Note first that if A^ = 2 the fixed point contributes nothing to the left side (from now 
on we drop the p subscript from N p ), so according to the formula we can have any number 
of fixed points of this type. These are called regular fixed points, and can only have a 
index of —1 (N/2 = 1, so the only permutation of sectors is the identity). Then, ignoring 
regular fixed points, the only possibilities for singularities are either one with A^ = 6, or two 
with N = 4 (recall A^ must be even). 

First consider the case of two singularities with A^ = 4. Since L(<p) = 1, at least one of 
the singularities must be fixed with its sectors permuted, since this is the only way to get 
a positive index. When A^ = 4 the permutation involves only two sectors, so <fi 2 will fix the 
sectors and the singularity will have an index of 1 — A^ = —3. However, the other singularity 
can contribute at most +1 to L(<f) 2 ), so this is incompatible with the Lefschetz number 
L((f> 2 ) = — 1. We conclude that this polynomial cannot correspond to a pseudo-Anosov 
mapping involving two singularities with A^ = 4. 
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Now consider the second case: one singularity with N = 6. Since L(<p) = 1, the sectors 
must be permuted to get a positive index. When N = 6 the permutation involves N/2 = 3 
sectors (Figure [5j. Therefore, 2 also permutes the sectors (so that the singularity con- 
tributes +1 to the Lefschetz number) and 3 fixes the sectors (so the singularity con- 
tributes 1 — N = — 5). Then, in order to match the Lefschetz number L(<p 2 ) = — 1 there 
must be two regular fixed points, which together make a regular period-2 orbit. Regular 
fixed points must appear at a multiple of their period (that is, a period-fc orbit becomes k 
fixed points of 4> k , (j) 2k , etc.). Hence, we conclude that a period-2 regular orbit exists, so that 
L((p 2 ) = 1 — 2 ■ (—1) = —1, as required. 

Next, for the third iterate 3 , the Lefschetz number is L((p 3 ) = —5: this is fine since the 
3 sectors of our N = Q singularity are now unpermuted, and the index is thus 1 — N = 
1 — 6 = —5. This means that there are no regular orbits of period 3. For the fourth iterate, 
we have £(0 4 ) = —5 again, but this time our main singularity has its sectors permuted, so 
its index is +1. Our earlier period-2 orbits appear at this iterate and contribute —2, and 
that leaves a deficit of —4, which is accounted for by invoking a single period-4 orbit. 

Let's look at one final iterate: for the fifth one L(<p 5 ) = —14, so something big is hap- 
pening. Our main singularity with N = 6 still has its sectors permuted, so its contribution 
is +1. No other periodic orbits that we've encountered before can contribute, since 5 is a 
prime iterate. Hence, the deficit of —14 — (+1) = —15 must be entirely made up by three 
period-5 orbits. 

So far the tally is (i) a single fixed point with N = 6; (ii) a regular (N = 2) period-2 
orbit; (iii) a regular period-4 orbit; and (iv) three regular period-5 orbits. Of course, we can 
keep going, but we have made our point: the list of Lefschetz numbers allows us to deduce 
which periodic orbits must appear at each iterate. Moreover, it is clear that the list is not 
arbitrary: there are many sequences of Lefschetz numbers that would be inconsistent, in the 
sense that no compatible periodic-orbit structure exists. What is remarkable is that all this 
information can be obtained from a simple polynomial. 

The littlest pseudo- Anosov 

The tools discussed so far allow us to answer an interesting question: for a given surface 
of genus g, what is the 'simplest' pseudo-Anosov mapping? By simplest, we mean the one 
with the least dilatation A, which we call 5+ (the '+' superscript corresponds to MCG + ). 
This is a question that has occupied some low-dimensional topologists for 20 years, ever 
since Penner first posed it. Originally it wasn't clear how this 'minimizer' should behave 
with the genus, but better and better examples were constructed, and currently the best 
general bound stands at <5+ ^ (2 + v^3) 1//s , implying that the minimizer converges to unity for 
large genusP^ In fact, this bound is sharp for g = 1^ and 2j n * 13 * 16 * ^ ^ but explicit pseudo- 
Anosov mappings have been found with dilatation below the bound for 3 < g < 5 J 1 1 * 1 9 * 
In genus 5, the pseudo-Anosov mapping associated with 5§ has the mysterious Lehmer's 
number ~ 1.17628 as a dilatationJ-^EQJ Moreover, the minimum dilatations for genus 2 to 5 
are all Salem numbers: these are roots of polynomials with only one eigenvalue outside the 
unit circle in the complex planed Nobody knows what happens in general for higher genus. 

In Ref. [ED we used a very direct method to prove, amongst others things, the minimality 
of the Salem numbers 6£ ~ 1.40127, 6% ~ 1.28064, and 6$ ~ 1.17628. Here, we illustrate the 
method for genus 3: Start with the polynomial associated with the candidate pseudo-Anosov 



6 



mapping, whose dilatation is believed to give 5% , in this case 

P(x) = x 6 - x 4 - x 3 - x 2 + 1 

with dilatation ~ 1.40127. Now how many degree 6, reciprocal, integer-coefficient polynomi- 
als have a largest real root that is less than this dilatation? A simple computer search tells us 
that the list is very small: there is only one, which can be factored as (x 3 — x — l)(x 3 + x 2 — l), 
with Lefschetz numbers 

L(0) = 3, L(0 2 ) = -1, L(4> 3 ) = -3, L(0 4 ) = 3, L(0 5 ) = -7, L(0 6 ) = -1 . . . 

where is the hypothetical pseudo-Anosov mapping. Since there are a finite number of 
possible singularities (because of the Euler-Poincare formula with % 3 = —4), we can check if 
this sequence of Lefschetz numbers is possible. Because L(<p) = 3, must fix at least three 
singularities. For instance, assume that fixes one singularity with N = 6, permuting its 
sectors as in Figure [5j and two with N = 4 with their N/2 sectors permuted. Then <p 2 will fix 
the sectors of the two singularities with N = 4, and each singularity will contribute 1—4 = —3 
to the Lefschetz number, for a total of —6. But from above, L(<f) 2 ) = —1, which is not 
possible: there is no way to make the Lefschetz number more positive by adding regular 
orbits, since they contribute negatively. Hence, <ft cannot fix these three types of singularities. 
There is only one other possible case — four singularities with N = 4 — which we can rule 
out as well. We conclude that the polynomial (x 3 — x — l)(a; 3 + x 2 — 1) is not associated 
with any pseudo-Anosov mapping. Thus, we have the mimimizer 5% ~ 1.40127. In order to 
truly complete the proof, one needs to explicitly construct the action of this pseudo-Anosov 
mapping on the surface (see Ref. ITT]) P^ 

In the end, is this numerological exercise worth it? Well, it certainly simplifies the study 
of pseudo-Anosov mappings on surfaces. In fact, even these abstract surfaces of genus g 
have practical significance: pseudo-Anosov mappings of disks, which can be associated to a 
stirring device such as that in Figure [2j can be lifted to surfaces of genus g by a 'branched 
cover.' For such applications, the biggest challenge is to develop practical tools (such as 
train track^) that can handle singularities with an odd number of sectors, which are math- 
ematically more challenging. 
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